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Abstract 



1-^ Attention is given to surface waves of shear-horizontal modes in 

^H piezoelectric crystals permitting the decoupling between an elastic in- 

plane Rayleigh wave and a piezoacoustic anti-plane Bleustein-Gulyaev 
wave. Specifically, the crystals possess 4 symmetry (inclusive of 42m, 
, 43m, and 23 classes) and the boundary is any plane containing the 

J> normal to a symmetry plane (rotated y-cuts about the Z axis). The 

^^ secular equation is obtained explicitly as a polynomial not only for the 

QQ metallized boundary condition but, in contrast to previous studies on 

\^ the subject, also for other types of boundary conditions. For the 

■^ metallized surface problem, the secular equation is a quadratic in the 

^^ squared wave speed; for the un-metallized surface problem, it is a 

sextic in the squared wave speed; for the thin conducting boundary 

r* problem, it is of degree 16 in the speed. The relevant root of the secular 

. ^H equation can be identified and the complete solution is then found 

/S (attenuation factors, field profiles, etc.). The influences of the cut 

j^ angle and of the conductance of the adjoining medium are illustrated 

numerically for GaAs (43m), BaLaGaaOy (42m) and Bii2Ge02o (23). 
Indications are given on how to apply the method to crystals with 222 
symmetry. 



1 Introduction 

The nature and properties of a piezoacoustic surface wave depend heavily 
on the crystallographic and anisotropic properties of the piezoelectric sub- 
strate, on the direction of propagation, and on the orientation of the cut 
(boundary) plane. For certain choices, the in-plane components of the me- 
chanical displacement decouple from the anti-plane component, leading to 
two different types of surface waves namely, the Rayleigh wave, elliptically 
polarized in the plane containing the direction of propagation and the nor- 
mal to the substrate surface, and the shear-horizontal (SH) wave, linearly 
polarized in the direction normal to the direction of propagation and parallel 
to the free surface. Moreover, either wave or both waves may be coupled 
to the electromagnetic fields. Of special interest are the configurations al- 
lowing for a piezoacoustic SH wave, decoupled from a purely elastic Rayleigh 
wave. Indeed, the former type of wave, also known as Bleustein-GulyaevPQI^] 
wave, penetrates more deeply into the substrate than the latter type; conse- 
quently, the acoustic energy is less localized and the power can be increased 
significantly before damage occurs (Tseng [3]) although, as pointed out by a 
Referee, there are exceptions to this behavior [3]. Moreover, a SH wave-based 
resonator is smaller than a resonator based on the propagation of a sagittally 
polarized surface wave (see Kadota and collaborators |S1IS]); this feature re- 
sults in a downsizing of design, an attractive point for the miniaturization of 
mobile phones for instance, where SAW devices are used as filters. Koerber 
and Vogel[71[8] identified all the cuts and rotations of axes leading to piezoa- 
coustic SH modes; they exist for some suitable cuts and transformations in 
the following crystal classes: 2, 23, 43m, 222, 2mm, 4, 4, 6, 4mm, 6mm, 
32, 42m, 6m2, 422, and 622. The main purpose of this paper is to derive 
explicitly the secular equation for piezoacoustic SH surface waves uncoupled 
from purely elastic Rayleigh waves polarized in a plane of symmetry, for a 
crystal in the 4 symmetry class (and thus for the 42m, 43m, and 23 classes). 
Several workers addressed this topic in the wake of the seminal papers 
by Bleustein and by Gulyaev, but explicit results remained limited either to 
propagation in special directions for which one of the piezoelectric constant is 
zero [31 m Uni [m [121 us] or to the case where the free surface of the substrate 
is metallized [9], or in the weak piezoelectric coupling approximation [9l [Ti] . 
Here the crystal is cut along a plane containing the Z axis and making any 
angle with the XY crystallographic plane. Moreover the surface may be 
metallized, or in contact with the vacuum, or in contact with a thin conduct- 



ing layer with arbitrary finite conductance. Also, no approximation is made 
about the strength of the piezoelectric effect. Attention is however limited to 
crystals with tetragonal 4 symmetry (inclusive of the tetragonal 42m, cubic 
43m, and cubic 23 symmetries). This limitation is not essential but simplifies 
the notation to a certain extent; it is nevertheless possible to extrapolate the 
method presented hereafter to crystals with lower symmetries such as or- 
thorhombic 222, as pointed out at the end of the paper. Note that the study 
of piezoacoustic SH surface waves in Potassium Niobate (2mm symmetry) 
was recently undertaken by Mozhaev and Weihnacht [I5j who showed that 
for that class, the secular equation is a cubic in the squared wave speed. 

The constitutive equations and the piezoacoustic equations for the class 
of crystals listed above are recalled in II. A and II. B, respectively, and some 
fundamental equations, which encapsulate the whole boundary value problem 
and its resolution, are quickly derived in II. C 

These fundamental equations are applied in Section III to the considera- 
tion of a piezoacoustic SH surface wave. A method based on the resolution 
of the propagation condition for partial modes first, and the resolution of 
the boundary value conditions next, would lead to quite an involved anal- 
ysis, including the analytical examination of a quartic polynomial for the 
coefficients of attenuation. The method based on the fundamental equations 
derived in II. C circumvents the stage of the quartic and delivers directly the 
secular equation in polynomial form. In particular it is seen that for metal- 
lized (a.k.a. short-circuit) boundary conditions (III. A), this equation is just 
a quadratic in the squared wave speed, whose relevant root is readily iden- 
tified. For open-circuit boundary conditions (III.B), the secular equation is 
also a quadratic in the squared wave speed, but it is not valid (the corre- 
sponding solution does not satisfy the boundary conditions.) For the free 
(non- metallized) boundary condition (III.C), the secular equation is a sextic 
in the squared wave speed. For the conductive thin layer boundary condition 
(III.D), the secular equation is a polynomial of degree 16 in the wave speed. 

Once the secular equation is solved for the speed, the complete description 
of the wave follows naturally (IV. A), including the attenuation coefficients 
and the profiles for the mechanical displacement, electric potential, traction, 
and electrical induction. A simple check for the validity of the solution is 
proposed in IV. B. 

The results are illustrated numerically and graphically using experimen- 
tal data available for GaAs (cubic 43m), BaLaGasOy (tetragonal 42m) and 
Bii2Ge02o (cubic 23) in V. The range of existence of the free SH wave with 




Figure 1: Rotated Y-cut about the Z axis. 

respect to the angle of cut, the speeds of propagation, the amphtude of the 
profiles, etc. are all quantities which can be obtained numerically with as 
high a degree of numerical accuracy as is needed. 

The paper concludes with a discussion of the merits and possible appli- 
cations of the method presented, and with a chronological account of the 
several advances in the field, without which the results of this paper could 
not have been established. 

2 Preliminaries 



2.1 Constitutive equations 

Consider a piezoelectric crystal with mass density p, possessing at most the 
tetragonal 4 symmetry (this symmetry includes the tetragonal 42m, cubic 
43m, and cubic 23 cases.) Let Cijki, iijk, and Qj be its respective elastic, 
piezoelectric, and dielectric constants with respect to the coordinate system 
XYZ along the crystallographic axes. Now cut the crystal by a plane con- 
taining the Z axis and making an angle 6 with the XY plane. The new 
coordinate system a;ia:2a:3 (say), obtained after rotation of XYZ about Z, is 
defined by 



(1) 



and so, the plane of cut is defined by X2 = 0, as shown on Figure 1. 

In the new coordinate system, under the electrostatic approximation for 





COS 6 


sin^ 




X 


= 


-sin^ 


cos^ 




Y 







1 




Z 



the electrical field, the stress tensor components 0",^ and the electric induction 
components Di are related to the gradients of the mechanical displacement 
u and of the electrical potential cf) by the constitutive relations, 



<^ij — CijklUi^k + ^ijk4>,k) Di — CiklUi^k — ^ik4>,ki 



(2) 



where the comma denotes partial differentiation (here with respect to the Xk 
coordinates). Using the Voigt contracted notation for the c, e, and e, these 
relations are written in matrix form as 
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(3)' 
Exphcitly, the Cij and ej, are deduced from the Cj, and Cj,- in XFZ by 
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well-known relationships. In particular. 



C44 — C44, en — en, 

ei4 = ei4 cos 29 — 615 sin 29, 

ei5 = ei5 cos 29 + 614 sin 29. 



(4) 



2.2 Piezoacoustic equations 

Now consider the propagation of an anti-plane (SH) inhomogeneous wave in 
the half-space X2 > 0, travelling with speed v and wave number k in the xi 
direction, with attenuation in the X2 direction. It is known [16] that for the 
crystals under consideration this wave decouples entirely from its in-plane 
counterpart, a purely elastic two-component Rayleigh wave. Thus the wave 
is modelled as Ui = U2 = 0, and 



{^3,0} = {Usikx2),v{kx2)}e''^''''-^'\ 



(5) 



for some yet unknown functions f/3 and ip of kx2. Accordingly, the consti- 
tutive equations Eq. pi) lead to similar forms for the stress and electrical 



components, 

{a,„ A} = ik{Uj{kx2),diikx2)}e"'^'''-^'\ (6) 

where tn = t22 = tu = d^ = 0, and 

^33 = ciaf/i - iciiU^, 

h2 = -«C44^3 + ei4V? + ^eisv?', 

t-si = CmUs + eisv? - ^ei4V?', 

c^i = eisf/s - ieut/g - env^, 

c^2 = euUs + zeisf/g + ienip' . (7) 

Here and hereafter, the prime denotes differentiation with respect to the 
variable kx2- 

Using the above- introduced functions f/3, (p, tij, di, the classical equations 
of piezoacoustics, 

o-ij,j = pui^tt, Di^i = 0, (8) 

reduce to 

- ti3 + ^t'm = -XUs, -rfi + z4 = 0, (9) 

where X := pv"^. 

When the region 0:2 < above the crystal is the vacuum (permeability: 
eo) and the surface 0:2 = remains free of tractions, then the boundary value 
problem corresponding to Eq. ^ is that of piezoacoustic (Bleustein-Gulyaev) 
SH surface waves. 

2.3 Fundamental equations for the resolution 

The method of resolution rests on the property that the piezoacoustic equa- 
tions ([9]) can be written in the form, 

^' = zN^, where ^(A;x2) = [f/3, <^, ^32, ^2]^, (10) 

and N is a real 4x4 matrix which, brought up to any positive or negative 
integer power n has the following block structure. 



N" 



n(") n^") 



with K(") = (K("))T, N^") = (N(")f, (11) 



K(") (Ni"))T 
where """ denotes the transpose, and where the submatrices of N" are 2 x 



2 matrices. The first-order differential form Eq. (10) of the piezoacoustic 



equations dates back to Kraut [H], and before that, to Strohp^H] and others 
(see Fahmy and Adler[19j for references) for the purely elastic case. The 
subsequent analysis below builds upon several crucial contributions, which 
are listed and put into context at the very end of this article. 

Now, because the wave amplitude must vanish away from X2 = 0, ^{kx2) 
is such that 

^(oo) = 0. (12) 

Clearly, pre-multiplication of N" by I defined as. 
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where 1 
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(13) 



produces a symmetric matrix. 
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(n) 



n\T 



(IN'^) 



(14) 



SO that taking the scalar product on both sides of Eq. (10) by IN"^, where 
the overbar denotes the complex conjugate, leads to $, ■ IN"^' = i$, ■ IN"+^^, 
the right hand-side of which is purely imaginary. Taking the real part and 
integrating yields $, ■ IN"^ 



const. = (by Eq. (12)), and in particular 



^(0) ■ IN"^(0) = 0. 



(15) 



These fundamental equations [201 [2T| 122] allow for a completely analytical 
derivation of the secular equation, for a great variety of boundary condi- 
tions. Note that because N is a 4 x 4 matrix, there are at most three inde- 
pendent fundamental equations according to the Cayley-Hamilton theorem. 
Any choice of three different integers n is legitimate, although the choice 
n = —1,1,2 seems to yield the most compact expressions for the compo- 
nents of N". 



3 Secular equations 

For the piezoacoustic (Bleustein-Gulyaev) shear-horizontal wave, the matrix 



N in Eq. ( 10 ) is written in compact form using the following quantity. 
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Note that N is indeed of the form Eq. (11). To express exphcitly integer 



powers of N, it proves convenient to use scalar multiples of the matrix IN", 



for which the fundamental equations Eq. (15) are also valid. For instance, 



with M(^\ M*- ■*■) respectively, defined by 



the fundamental equations (|15|) also hold at n = 2, — 1 when IN" is replaced 
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These symmetric matrices are given explicitly by their components 
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and 
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m(,-^) = (1 + '-^)^'X - C44(l + n\ m[-'^ = -'-^^'X + 
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"11 --14 

Now all the required equations and quantities are in place to treat various 
electrical boundary value problems. 

3.1 Metallized (short-circuit) boundary condition 

Here the surface of the crystal is coated with a thin metallic film, with thick- 
ness negligible when compared to the wavelength, and brought to a zero 
electrical potential. Moreover, the coating still allows the surface to remain 
free of mechanical tractions. Then 



(J23 = 0, = 0, at Z2 = 0, so that ^(0) = [UsiO), 0, 0, c?2 

Writing ^(0) = f/3(0)[l, 0, 0, a], where a = d2iO)/Us{0) is complex, the fun 
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lead to the following homogeneous system of equations. 
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(22) 



For a non-trivial solution to exist, the determinant of the system's matrix 
must be zero. Factoring out common factors, this condition reads 



X — 044(1 + ft^) C44K^ — C44 

2[X-2c44(l + /«^)] -X + 2c44(l + 2fi;2) -4c44 

;i + ^-^)k^X - C44(l + K^) {X - C44)fi:2 X - C44 



'14 



0. (23) 



When X = 0, the first column in the determinant becomes proportional to 
the third, and the determinant is zero. Hence X is a factor of the determi- 
nant; the remaining factor is a quadratic in X, 

X- - C44( ''"'" + '^^^ + ^'h x + 24( ^"'" + ''^^ t ^"h = 0, (24) 
eiiC44 + efg eiiC44 + efg 

that is, the explicit secular equation for piezoacoustic (Bleustein-Gulyaev) 
anti-plane surface waves on a metallized tetragonal 4 (or tetragonal 42m, or 
cubic 43m, 23j crystal, cut along any plane containing the Z axis. 

This equation being a quadratic in X, it is solved explicitly and it yields 
a priori two roots. The selection is made by considering the known speed of 
a Bleustein-Gulyaev surface wave in a cubic 43m or 23 crystal when 6 = 45°. 
Then the root of the quadratic corresponding to the plus sign is X = pv'^ = 
2c44, and the root corresponding to the minus sign is 

2 ^2 

PV' = C44(l + ^^^) = C44(l + ^^^), (25) 

044611 + e^^ £44611 + Ci^ 

in accordance with Tseng ^ (see also Koerber and Vogel|7j, Alburque and 
Chao[TO], Velasco[TT], Bright and Hunt|T3]). By continuity with the other 

cases (4, 42m, 43m, 23, 6 7^ 45°), the speed of the generic Bleustein-Gulyaev 



wave found from the secular equation (24) is fBCm, given by 

P^Gm ^ 3611C44 + 4e^4 + 4e^5 - [(eiiC44 + ^ej^f + (4614615)^)2 
C44 2(611644 -|- 615) 

For tetragonal 42m, cubic 43m, or cubic 23 crystals, 615 = and by Eq. (|4]), 
this expression reduces to 

P^BGm 3 + 4x^ - [(1 + Ax^ cos^ 2^)2 + 4x^ sin^ 46] \ 



644 2(l + x2sin2 20) 



(27) 



as proved by Braginskii and Gilinskii[n] using a different method. In Eq. (27), 



X^ = 614/(611644) is the "piezoelectric coupling coefficient" for bulk waves. 

3.2 Electrically open boundary condition 

The substrate is said to be "mechanically free, electrically open" (Ingebrigsten 
Lothe and Barnett|24l) when 



^32 = 0, D2 = 0, at X2 = 0, so that ^(0) = [[73(0), (/?(0), 0, 0]^. (28) 

10 



Writing ^(0) = f/3(0)[l, a,0, 0], where a = (/?(0)/f/3(0) is complex, the 
fundamental equations (15) for IN {n = 1), M*^^-' (n = 2), and M*^^-'-) (n = 
— 1) lead to the following homogeneous system of equations, 
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(29) 



For a non-trivial solution to exist, the determinant of the system's matrix 



must be zero. Its components are given in Eq. (17), Eq. (19), and Eq. (20). 
As in the short-circuit configuration, at X = the first and third columns in 
the determinant become both proportional to [1,4, 1]""" and so X is a factor 
of the determinant; the remaining factor is a quadratic in X, 
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(30) 

At this stage, an important point must be raised. Although this secu- 
lar equation might seem legitimate at first sight, it must be recalled that 



it was obtained through a process based on the fundamental equations ( 15 ) 



which, due to the involvement of integer powers of the matrix N, might gen- 
erate spurious secular equations. In fact, it has been proved [211 [H] that the 
Bleustein-Gulyaev wave does not exist for open circuit boundary conditions. 



Hence the secular equation Eq. (30) is not valid. It is given here for complete- 



ness and to illustrate one limitation of the fundamental equations approach. 
However, as is seen in IV.B, a simple check can be done to realize whether 
the speed given by an explicit secular equation is valid or not. 



3.3 Free boundary condition 

In the general case of a non-metallized, mechanically free boundary, the tan- 
gential component of the electric field and the normal component of the elec- 
tric induction are continuous across the substrate/vacuum interface. These 
continuities lead to the relationship (e.g. Dieulesaint and Royer [23] p. 288), 



^2(0) = ieo¥?(0), so that ^(0) = ip{0)[a,l,0,iec 



(31) 



where a = ai+ia2 = U3{0)/ip{0) is complex. Now the fundamental equations 
([T5I) for IN (n = 1), M(2) (n = 2), and M^-^^ {n = -1) lead to a non- 
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homogeneous linear system of equations, 

N32 eoN2i N31 
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By Cramer's rule, the unique solution to this system is 

2ai = Ai/A, 2a2 = A2/A, a? + a^ = A3/A, 



(32) 



(33) 



where A is the determinant of the 3x3 matrix in Eq. (32), with components 
given in Eq. (17), Eq. (19), and Eq. (20), and the A^ are the determinants 



obtained by replacing this matrix's A;-th column with the vector on the right 



hand-side of Eq. (32). It follows from Eq. (33) that 

A? + A2-4AA3 = 0, 



(34) 



which is the explicit secular equation for piezoacoustic (Bleustein-Gulyaev) 
anti-plane surface waves on a non-metallized, mechanically free 4 (or 42m, 
4:3m, 23^ crystal, cut along any plane containing the Z axis. 

The expansions of the determinants A, . . . , A3 are lengthy and are 
not displayed here, but they are easily computed by using the components 



Eq. (17), Eq. (19), and Eq. (20). It turns out that A factorizes into the prod- 
uct of eo and a cubic in X which is independent of eo, while Ai (resp. A2, 
A3) factorizes into the product of CqX (resp. X, eo) and a polynomial which 



is quadratic in X and linear in e^. The resulting secular equation (34) is a 



sextic in X and a cubic in eg, with the coefficient of the e^ term proportional 



to the "metallized secular equation" (24) and the coefficient of the ef] term 



proportional to the (non- valid) "open circuit secular equation" (30). 

It is emphasized again that, as in III.B, great care must be taken to 



ensure that the speed given by the secular equation Eq. (34) leads to a valid 
solution. This point is discussed in IV. B. 



3.4 Thin conducting layer boundary condition 

As a final type of boundary condition, consider that the semi-infinite sub- 
strate is covered with a metallic ffim with thickness h and conductance 7, 
where h is assumed to be so small with respect to the acoustic wavelength 
that the effects of mechanical loading can be neglected. Then the perme- 
ability of the region X2 < close to the interface is changed from eo (see 
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previous Subsection) to cq — ie (see Royer and Dieulesaint[2S], P- 301), with 
e = 'jh/v. Replacing the former quantity with the latter in the previous 
Subsection leads to similar results as above, except that the 3x3 matrix and 



the right hand-side of Eq. (32) are now replaced with 



N32 + eN2i 
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1) _ 


26Mt^) 


- (4 + e^)Mt,'^ 



respectively. Then the secular equation is Eq. (34) where A 



, (35) 



. , A3 are 

appropriately changed. Because here the quantity e depends upon v, the 
secular equation is a polynomial in the speed of degree higher than in the 
previous Subsection, namely it is a polynomial of degree 16 in v. 

4 Construction of the solutions 

4.1 Description of the wave 

Once a wave speed is determined from the secular equation, it is a rather 
straightforward matter to construct the corresponding complete solution. 
Indeed, the anti-plane mechanical displacement M3, the electrical potential 
0, the shear stress (T32, and the electric induction D2 are given by 

[W3, 0,^32, 1^2] (xi,X2,t) = ^{[U3,ip,^kts2,^kd2]{kx2)e''^''''^'^}. (36) 

Here U3, ip, ^32, ^2 are the components of $,, solution to $,' = iN^. Taking ^ 
in exponential form leads to the following decaying solution, 

^ = /Ji^ie'^^"!"^ + 132^6''"''^' , (37) 

where /3i, (32 are constants, gi, g2 are the two roots with positive imaginary 
part to the inhomogeneous wave propagation condition: det (N — ql) = 0, 
and the ^ satisfy: N^* = gj^*. 

Explicitly, the qi are roots of the quartic, 

(eiiC44 + 615)^"^ - 4ei4ei5g^ - [eii(X - 2C44) - 4ei4 + 2el^]q^ 

+ 4ei4ei5g + ei5-eii(X-C44) = 0, (38) 
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and the ^* are proportional to any column vector of the matrix adjoint to 
N — gjl, the third one say. Hence, 

-\ T 



e 



C44 C44 



(39) 



where the non-dimensional real quantities ai, bi, fi, gi {i = 1,2) are given by 
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(40) 



Finally, the ratio /32//3i comes from the condition that the third component 
of ^(0), proportional to ^23(0), must be zero, so that /32//3i = — /i//2- 



4.2 Validity of the solution 

Once a wave solution has been constructed, its validity must be checked that 
is, it must satisfy the boundary conditions. 

Thus, for a solution to the short-circuit problem (III. A), it must be 
checked that ip{0) = f3ibi + /32b2 and ^32(0) = /3iei4/i + (3261^/2 are indeed 
equal to zero. These conditions are equivalent to checking that 
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(41) 



For a solution to the open-circuit problem (III.B), it must be checked 
that t32(0) = /3iei4/i + /32ei4/2 and ^2(0) = (3iiel^/cu)gi + P2{elJcii)g2 are 

/l /2 



equal to zero. These conditions are equivalent to checking that 
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0. 



As expected [Ml E], this condition is never satisfied. 

For a solution to the free boundary problem (III.C), it must be checked 
that ^32(0) = PiCiifi + (32eiif2 = and that ^2(0) = zeoV5(0) that is, 
/3i (614/044)5^1 + /92(ei4/c44)5'2 = i^.Q^Pibi + /92&2)- These conditions are equiv- 
alent to checking that 
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(42) 
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In general, this condition is met only for a limited range of the cut angle 9. 
Finally, for a solution to the thin conducting layer boundary problem 
(III.D), it must be checked that ^32(0) = and that c?2(0) = (e + zeo)<^(0). 
These conditions are equivalent to checking that 



(43) 
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5 Examples 

In this Section, the secular equations derived in III and the tests presented 
in IV are employed to find numerically the wave speed and its range of exis- 
tence for three crystals, one with cubic 43m symmetry, one with tetragonal 
42m symmetry, and one with cubic 23 symmetry. Data collected from the 
specialized literature are used for the values of the mass densities, of the stiff- 
nesses, and of the piezoelectric and dielectric constants. In order to graph 
the depth profiles, a frequency of 100 MHz and a mechanical displacement 
of 10"^'^ m at a;2 = are picked, to fix the ideas. Note that at a 45° angle of 
cut, the profiles present pure (non-oscillating) exponential decay, because the 



propagation condition (38) is a biquadratic and the corresponding roots are 
purely imaginary; they are essentially similar to those displayed by Bright 
and Hunt[13j. Here the profiles are computed at angles 7^ 45°. 

5.1 AlAs 

For Aluminum Arsenide (43m symmetry) the physical quantities of interest 
arei27j: p = 3760 kg-m'^, C44 = 58.9 x 10^ N-m'^, 614 = -0.225 C-m'^, and 
en = 10.06eo. 

Using the results of HI and IV, it is found that the speeds of the piezoa- 
coustic SH surface wave with metallized (HI. A) and with free (III.C) bound- 
ary conditions are almost indistinguishable on a graph from the speed of the 
bulk shear wave. For instance at 45°, these speeds are (m-s~^): 3976.784, 
3976.965, and 3976.966, respectively. Note however that the SH surface wave 
for the metallized ("shorted") boundary condition exists for all values of 6 
(within the range delimited above by the speed at 45° and below by the 
speed at 0° and at 90°, which is 3957.890 m-s~^) whereas the SH surface 
wave for the un-metallized ( "free" ) boundary condition exists only within a 
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Figure 2: Speeds of piezoacoustic waves in a AlAs (43m symmetry) crystal, 
as a function of the cut angle: bulk shear wave (upper curve), SH surface 
wave for free (un-metallized) boundary conditions (intermediate curve), and 
SH surface wave for metallized boundary conditions (lower curve). 

limited range, delimited above by the speed at 45° and below by the speed 
at 45.0° ± 8.51°, which is 3975.276 m-s'^ 

Figure 2 displays the variations of the three speeds as a function of 9. The 
speed of the bulk shear wave is always above the speed of the SH surface wave 
for the metallized boundary condition; they are both defined everywhere. 
The speed of the SH surface wave for the un-metallized boundary condi- 
tion is intermediate between these two speeds, but exists only in the range 
[36.49°, 53.51°]. A zoom is provided for this range. In that zoom, the curve 
for the bulk shear wave almost coincides with the curve for the SH surface 
wave corresponding to the un-metallized boundary condition; together they 
form the upper curve whilst the lower curve represents the variations of the 
SH surface wave speed corresponding to the metallized boundary condition. 

Note that the simple test for the solution's validity presented in IV.B 
works perfectly here and henceforward. Thus using a 40 digit precision under 



MAPLE for AlAs, the modulus of the determinant in (42) is found to be less 
than 10-19 at 45.0° ± 8.5125° and more than 0.2 at 45.0° ± 8.5129°. 

Figure 3 shows the variations with depth of the fields of interest (me- 
chanical displacement, shear stress, electrical potential, electric induction) 
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Figure 3: Depth profiles of tlie SH surface wave for metallized boundary 
conditions in a AlAs crystal cut at 22°30': mechanical displacement, shear 
stress, electrical potential, electric induction. 

for the SH surface wave corresponding to the metallized boundary condition 
at ^ = 22°30'. The variations of the fields are presented over 250 wavelengths, 
and zooms are provided for the [0, 5] wavelengths range, where 0, (J32, and 
D2 undergo rapid changes. 

5.2 BLGO 



Soluch et a/.[2H] measured experimentally the elastic, piezoelectric, and di- 
electric properties of BaLaGa307 (42m symmetry, mass density: p = 5450 
kg-m-3) as: C44 = 39 x 10^ N-m'^, eu = 0.29 C-m'^, and en = 12.4eo. 

Here, the speeds of the piezoacoustic SH surface wave with metallized 
(HI.A) and with free (III.C) boundary conditions differ more notably than in 
the previous example from the speed of the bulk shear wave. For instance at 
45°, these speeds are (m-s-^): 2700.739, 2701.239, and 2701.242, respectively. 
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Figure 4: Speeds of piezoacoustic waves in a BaLaGasOr (42m symmetry) 
crystal, as a function of the cut angle: bulk shear wave (upper curve), SH sur- 
face wave for free (un-metallized) boundary conditions (intermediate curve), 
and SH surface wave for metallized boundary conditions (lower curve). 



The range of values for the wave speed fBCm of the metallized "shorted" 
boundary condition is delimited above by the speed at 45° and below by 
the speed at 0° and at 90°, which is 2675.063 m-s~^. For the un-metallized 
"free" boundary condition, the corresponding (limited) range for the wave 
speed fBGf is bounded above by the speed at 45° and below by the speed 
at 45.0° ± 7.485°, which is 2699.369 m-s~^. The difference between the two 
speeds is the largest at 6^ = 45°; there the ratio 2(t'BGf — 'i^BGm)/'yBGf is equal 
to 3.70 xlO"! 

Figure 4 displays the variations with 9 of the speeds for the bulk shear 
wave, for the SH surface wave corresponding to the un-metallized bound- 
ary condition, and for the SH surface wave corresponding to the metallized 
boundary condition. Figure 5 shows the variations with depth of the fields of 
interest (mechanical displacement, shear stress, electrical potential, electric 
induction) for the SH surface wave corresponding to the metallized boundary 
condition dX 9 = 22°30'. Similar comments to those made for Figure 3 apply. 
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Figure 5: Depth profiles of tlie piezoacoustic SH surface wave for metal- 
lized boundary conditions in a BaLaGasOy crystal cut at 22°30': mechanical 
displacement, shear stress, electrical potential, electric induction 

5.3 BGO 

The relevant physical quantities of Bismuth Germanium Oxide (Bii2Ge02o, 
23 symmetry) arefSH]: P = 9200 kg-m'^, C44 = 25.52 x 10^ N-m'^, e^ = 0.983 
C-m~^, and en = 38.0eo. 

Here the differences between the speeds of the bulk shear wave, of the 
SH surface wave corresponding to the metallized boundary condition, and of 
the SH surface wave corresponding to the un-metallized boundary condition 
are more marked than in the previous example. At 45°, these speeds are 
(m-s-i): 1747.812, 1756.836, and 1756.846, respectively. At 0° and at 90°, 
the speeds of the bulk shear wave and of the SH surface wave corresponding 
to the metallized boundary condition are both equal to 1665.507 m-s~^. The 
SH surface wave for the un-metallized boundary condition exists only in the 
range 45.00±3.149°, and at the extremities of this range, its speed is 1755.068 
m-s~^. 

Figure 6 displays the variations of the three wave speeds as a function of 
9. Figure 7 shows the variations of the quantity 2(t;BGf — 'yBGm)/'yBGf with the 
angle of cut; its largest (smallest) value is 1.027x10"^ at 45.00° (0.975x10-^ 
at 48.149°). 
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Figure 6: Speeds of piezoacoustic waves in a Bii2Ge02o (23 symmetry) crys- 
tal, as a function of the cut angle: bulk shear wave (upper curve), SH surface 
wave for free (un-metallized) boundary conditions (intermediate curve), and 
SH surface wave for metallized boundary conditions (lower curve). 
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Figure 7: Variations of 2(t;BGf ^ 'i'BGm)/'yBGf with ^ in a Bii2Ge02o (23 sym- 
metry) crystal. 
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Figure 8: Depth profiles of tlie piezoacoustic SH surface wave for metal- 
lized boundary conditions in a Bii2Ge02o crystal cut at 22°30': mechanical 
displacement, shear stress, electrical potential, electric induction. 

Figure 8 shows the variations with depth of the fields of interest (me- 
chanical displacement, shear stress, electrical potential, electric induction) 
for the SH surface wave corresponding to the metallized boundary condition 
at ^ = 22°30'. Similar comments to those made for Figures 3 and 4 apply. 

Figure 9 shows the variations with depth of the same fields for the SH 
surface wave corresponding to the un-metallized boundary condition at ^ = 
42°30'. By comparison with the previous Figure it can be seen that the SH 
surface wave for the un-metallized boundary condition penetrates far more 
deeply than the SH surface wave for the metallized boundary condition. The 
electrical potential and the electric induction are plotted inside the crystal 
for the range [0, 500] wavelengths, and also in the vacuum over the crystal 
for the range [-10, 0] wavelengths; the continuity of these fields across the 
interface is made apparent with a zoom for the range [-5, 5] wavelengths. 



6 Concluding remarks 

The method of resolution for the title problem of this paper is based on the 



fundamental equations (15). This method has proved itself to be very effec- 
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Figure 9: Depth profiles of tlie piezoacoustic SH surface wave for free (un- 
metallized) boundary conditions in a Bii2Ge02o crystal cut at 42°30': me- 
chanical displacement, shear stress, electrical potential, electric induction. 
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tive and versatile. The end result is the complete analytical elucidation of the 
problem, for a great variety of surface impedance problems in a piezoelectric 
half-space i.e., problems where the electric induction is proportional to the 
electrical potential: k^ = ivZD2 (say) at the boundary plane. The method 
can be followed through when the impedance Z is zero (short-circuit), infinite 
(open-circuit), pure imaginary (free boundary), or complex (thin conducting 
layer) . It has already been used for other types of surface impedance prob- 
lems for elastic interface waves (Stoneley waves [20| 122] . Scholte waves [21]) 
and could be adapted to configurations [30j with a resistance force propor- 
tional to the normal velocity, a mass concentrated in a thin surface layer, 
a system of elastic oscillators resting on an elastic half-space, a thin elastic 
layer longitudinally deformable, etc. The method can also accommodate a 
coupling between elastic and piezoelectric fields, in situations such as the one 
treated here or for instance, the case of interface acoustic waves at a domain 
boundary |3Tj. 

Here, attention was restricted to Bleustein-Gulyaev waves in tetragonal 
4 piezoelectric crystals. The extension to the classes of orthorhombic 222 or 
monoclinic 2 crystals is straightforward and only requires the computation of 



the elements of the matrix N in Eq. (10). Their general expression is given 



for instance by Abbudi and Barnett[32]. As an illustration, they are now 
presented for rhombic 222 crystals. 

For such a crystal, the relevant non-zero piezoacoustic constants in the 
crystallographic coordinate system are C44, C55, 614, 625, en, and €22- In the 
coordinate system obtained after the rotation Eq. ([I]), they are 

C44 = C44 cos^ 9 + C55 sin^ 9, C45 = (£44 — C55) cos 9 sin 9, 

C55 = C55 cos^ 9 + C44 sin^ 9, ei2 = (£22 - en) cos 9 sin 9, 

en = en cos^ 9 + €22 sin^ 9, e22 = €22 cos^ 9 + en sin^ 9, 

ei4 = ei4 cos^ 9 — 625 sin^ 9, eis = (ei4 + 625) cos 9 sin 9, 

625 = 625 COS^ 9 - 614 siu^ 9, 624 = -615. (44) 



The equations of motion can be cast in the form Eq. (10) where the matrix 



N is defined by its 2 x 2 blocks N^^^ N^^\ K^^) in Eq. lK^. The components 
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of — N-j^ and N2 are given here by 
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Here, the quantity k^ is defined by 
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Finally, in guise of a Conclusion, the main relevant advances toward the 
full resolution of the problem presented in the paper are recapitulated. In 
the purely elastic case, the secular equation for Rayleigh waves polarized 
in a plane of symmetry was derived by Currie[33] who, using an algebraic 
approach based on the Stroh formalism, obtained the equations. 



U(0) ■ k(")u(o) = 0, 



(48) 



where U(0) is the mechanical displacement on the free surface. Although 
these equations are also valid in generally anisotropic crystals, his derivation 
of the secular equation for triclinic (no symmetry) crystals apparently leads to 
a trivial identity. This problem was later corrected by Taylor and Currie[3i] 
and by Taziev[35] (see also Ting[36])- In contrast to these approaches based 
on the formulation of the equations of motion as a first-order differential 
system for the displacement-traction vector, Mozhaev[37j wrote the equations 
of motion as a second-order differential system for the displacement vector. 



aV" - i/5U' + 7U = 0, 



(49) 
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where a, /3, 7, are real symmetric matrices. Then, using first integrals, he 
quickly derived the secular equation for orthorhombic crystals. Destrade|38j 
rewrote the equations of motion, this time in the form 

at" - 0t + 7t = 0, (50) 

for the tractions, where a, /3, 7, are real symmetric matrices. Adapting 
Mozhaev's first integrals, he re-derived (unaware of Currie's result) the sec- 
ular equation for Rayleigh waves polarized in a symmetry plane. He also 
mentioned (and the proof was later given in the review article by Ting|39j) 
that the method of first integrals could not be used for arbitrary anisotropy 



when the equations of motion are written as Eq. (49 ) or Eq. (50 ). Recently [20l 
|2T| [22] he made the connection between Currie's and Taziev's use of integer 
powers of the Stroh matrix N and Mozhaev's first integrals, as shown also 
here in II. C Note that Mozhaev and Weihnacht|15j were able to solve the 
problem of SH surface modes of a 2mm crystal using first integrals of the 
piezoacoustic equations written as a second-order differential system. 
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